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Abstract 

X-ray synchrotron emission tells us of the highest energy reached by accelerated 
electrons. In a few supernova remnants (SN 1006, G347.3—0.5) this is the dominant 
form of X-ray radiation, but in most it is superposed to the dominant thermal 
emission. Thanks to the spectro-imaging capability of Chandra and XMM-Newton, 
X-ray synchrotron emission has now been unambiguously detected in most young 
supernova remnants (Cas A, Tycho, Kepler). It arises in a very thin shell (a few 
arcsecs) at the blast wave. The thinness of that shell (much broader in the radio 
domain) implies that the high energy electrons cool down very fast behind the 
shock. The magnetic field that one deduces from that constraint is more than 100 
fj ,G behind the shock. 

Key words: Acceleration of particles, Magnetic fields, Cosmic rays, ISM: 
supernova remnants, X-rays 
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1 Introduction 


The current paradigm is that the bulk of the cosmic-rays (up to the “knee” at 
3 x 10 15 eV) are ac celerated at the blast waves g enerated by supernova explo¬ 
sions in our galaxy (Blandford and Eichler. 11987 ). Indeed supernova remnants 
(SNRs) are all non-thermal radio emitters, attesting of the presence of acceler¬ 
ated electrons at energies of 1 GeV or so in larger amounts than in the average 
interstellar medium. In most SNRs, the radio emission is limb-brightened, con¬ 
firming that those accelerated electrons originate at the shock rather than in 
a central pulsar. 
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The X-ray emission in those shell-type SNRs is usually mostly thermal (dom¬ 
inated by strong lines of heavy elements), but the synchrotron emission may 
extend up to the X- rays and contr i bute t o the emission as well. In several 
like SN 1006 ( Kovama et all il995 ) and G347.3—0.5 ( Kovama et al 


cases, 


19971), the synchrotron component actually dominates the X-ray emission. 


The slope of the X-ray synchrotron emission is always much steeper (a ~ 1.5) 
than that of the radio emission (a ~ 0.5). This indicates that the X-rays 
(emitted by electrons at energies > 1 TeV) probe the cut-off of the electron 
distribution, and can tell something about the limits of the shock acceleration 
mechanism. The open questions on which observations of X-ray synchrotron 
in SNRs may shed light are: 

• How efficient is cosmic-ray acceleration in SNRs ? What is the energy density 
of accelerated particles ? 

• What is the maximum energy of accelerated particles ? 

• How large is the magnetic field ? Is it very turbulent ? Is it amplified ? 


Unfortunately, not much is known about the synchrotron emission in the seven 
decades gap between the radio and the X-rays, partly because it is masked by 
other components (like infra-red emission by heated dust) and partly because 
it is simply very difficult to detect faint extended emission. Because of that 
ignorance, it is often assumed that the spectrum extends as a power law up 
to the c ut-o ff, even though non linear acceleration models predict some con¬ 
cavity ( Baring et al. . 19991). The synchrotron emission has been detected in 
Cas A in the near infra-red ( Jones et al. . 2003), but this is associated with the 
bulk of the radio emission (near the interface with the ejecta), not with the 
relatively faint synchrotron emission from behind the blast wave discussed in 
the following sections. 


In this review I report on recent observations of X-ray synchrotron emission 
in three young SNRs dominated by thermal emission from the ejecta (Cas A, 
Kepler and Tycho), and two older SNRs dominated by non-thermal emission 
(SN 1006 and G347.3—0.5). The inference is that the magnetic field just behind 
the blast wave is quite large (up to 200 pG). 


2 Observations at high angular resolution 


Before it could be spatially resolved, the X-ray emission behind the blast wave 
in young SNRs was expected to be the thermal emission of the shocked am¬ 
bient gas, distributed more or less evenly between the blast wave and the 
outer boundary of the ejecta. But the Chandra results challenged that pre¬ 
conception. Th e imag es of the con tinuum emiss i on (4 to 6 keV) in Cas A 


( Gotthelf et al. . 2001 ) and Tycho (. Hwang et ah . 2002 ) clearly show a very 
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thin X-ray rim at the blast wave. The apparent width of the rim is less than 
4" or 2 10 1 ' cm in Cas A, and appears to be twice lower at places. That X-ray 
rim runs all around both remnants, albeit in a broken manner (particularly in 
Cas A). 

The 4 to 6 keV continuum emission in Kepler shows a very similar structure. 
Sharp (although not nearly as sharp) rim s have also been observed in SN 


rot nearly as snarp) rim s nave also been observed m al 
20031 Bamba et all 20031) and G347.3—0.5 ( Lazendic et al 


1006 ((Long et al 

2Q03). In those two SNRs, older and dominated by non-thermal emission, and 


contrary to the younger three mentioned above, there is no clear boundary 
behind which the thermal emission from the ejecta dominates. In G347.3—0.5, 
the thermal emission is not even detected at all. Also, because the non-thermal 
emission dominates, the energy range where its morphology can be studied is 
not restricted to the 4 to 6 keV band (SN 1006 is very faint in that band). 


In Tycho, the rim clearly marks the outer boundary of the X-ray emission, so 
it is natural to think of it as a thin sheet covering the entire sphere, rather 
than a linear filament. This is supported by the observation of SN 1006, where 
the rim is definitely resolved, and much sharper outwards than inwards. The 
visual appearance of the rims in Cas A is much less regular. This is not a 
strong argument against the sheet geometry, however, because it is actually 
very reminiscent of the optical emission in older remnants like the Cygnus 
Lo op, w hich is generally interpreted as due to a wrinkled sheet of emission 
(Hester. 119871) . In keeping with this evidence, I will assume the rims are places 
where the sheet of emission at the blast wave is observed tangentially. 


As noted by Berezhko and Volk (2004J), since this is observed in projection 
the scale height of the spherical layer must be even smaller than the observed 
width. I argue in App. A that it should be 4.6 times smaller. This conclusion 
does not depend on the curvature radius of the sheet (it does not have to be 
equal to the SNR radius, it could even be negative - outwards), therefore it 
is quite robust. Exceptional geometrical conditions (like model B of Hester, 
19871) . normally associate a narrower rim with a larger brightness contrast 
between the rim and its surroundings. This should be testable statistically. 


The brightness contrast of the X-ray continuum between at the rim and behind 
the rim is very large (typically 5). It is consiste nt with the thin s heet model 
in Cas A ( Berezhko and Volk . 2004 ) and Tycho ( Volk et al. . 20051) . but leaves 
little room for anything else. This implies that most of the volume between 
the blast wave and the interface with the ejecta is actually X-ray dark (in the 
continuum emission). 


Another important observation is the nearly featureless nature of the spectrum 
(very faint lines are observed). In a thermal framework, this can be explained 


if ionization is very far out of equilibrium ( Hwang et ah . 20021 . on Tycho). But 
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of course it is very natural if the spectrum of the rims is non-thermal. It is 
important to note in that respect that the line emission is not peaked at the 
rims. The brightness in the lines is essentially the same in a region behind 
the rims as on the rims in Cas A and Tycho. This suggests that the faint line 
emission comes from the entire volume between the blast wave and the bright 
ejecta, contrary to most of the continuum emission. 


3 The nature of the X-ray emission behind the blast wave 


The observed geometry is inconsistent with thermal models in a uniform 
medium which predict emission everywhere up to the interface, with only 
a slight maximum at the blast wave. This region is full of hot gas, which can¬ 
not cool down efficiently at those low densities, so it has to be entirely X-ray 
bright. Such a sharp decline of the X-ray emission behind the blast wave could 
be due to a recent density increase of the ambient gas, but it would have to 
be an extraordinary coincidence that this happens exactly at the same time 
all around the remnants in both Cas A and Tycho. In Cas A (a type II SN), 
having just hit a spherical wind shell could create such a condition, but it 
would result in a much s lower shock, contrary to the observed proper motions 
( DeLanev and Rudni ck. 20Q.‘lil . 


Other known transient phenomena related to the thermal gas cannot explain 
the observations at all: 

• Heavy element ionization affects the lines rather than the continuum. 

• Electron heating by the ions predicts a hardening of the spectrum, but not 
a decrease in intensity. 

• Dust destruction results in increasing emission. 


The only other possible source of X-ray radiation is the accelerated particles. 
Non-thermal bremsstrahlung (by suprathermal electrons at 10 keV or so) could 
be an option. The density of targets (the thermal gas) does not decline steeply 
behind the shock, but collisional losses with the thermal electrons could be 
strong enough to get rid of the particles themselves as they are advected down¬ 
stream. A gross estimate of the energy loss time for suprathermal electrons at 
energy E^ e y in gas of electronic density is trv.ni = 4 .9 x 10 7 s E^y/n e . A 
more accurate formula may be found in Sturner et ah (1997). Around 5 keV, 
this could be consistent with the observed width in Cas A (corresponding to 
an age of 2.5 xlO 8 s rj 4, where r is the total compression ratio) for a down¬ 
stream gas density of 1.8 cnr 3 (r/4) _1 . However, non-thermal bremsstrahlung 
must be associated with brighter thermal emission. In Cas A, Tycho or Kepler, 
no such bright thermal emission is observed associated with the X-ray rims. 
With an interstellar column density X H of several 10 21 cm -2 , it is conceivable 
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that the thermal gas could hide at low temperature around 1 or 2 xlO 6 K. 
Such a low tem peratu re is actu a lly pr edicted by strongly non linear accelera¬ 
tion models ( Decourchelle et al. . 20001) . SN 1006, on the other hand, has much 
lower IVh — 7 x 10 20 cm -2 (jPubner et al. . 2002 ). The soft thermal emission 
(in particular the O K lines) is clearly seen there at the same place as the 
harder X-ray rims. However it is not brighter than the non-thermal emission 
at all, which eliminates the possibility that the non-thermal continuum could 
be a tail of the thermal emission. 

The most natural remaining option is synchrotron emission by high energy 
electrons. Here again, advection of the particles and the magnetic held (with 
only slight adiabatic losses) cannot explain the very sharp drop behind the 
shock. On the other hand, the particles may lose their energy radiatively fast 
enough as they are advected (so that their synchrotron emission is shifted 
below the X-ray range) t o expla i n a very thin emission region if the mag¬ 
netic held is large enough f Ballet . 120031: Vink and Lami ng. :2003h . In Sect. 4 I 
develop the consequences of that interpretation. 


The synchrotron interpretation implies also that the density (or the tempera¬ 
ture) of the ambient gas must be low enough that it does not contribute sig- 
nihcantl y to the X-ray emission. This is possible if n (shocked) < 0.15 cm 


-3 : 


m 


Kepler (jCassam-Chenai et all l2004al) . The constraint for Tycho is probably 


a little less severe. Another prediction is that the radio emission should be 
much broader (not so peaked at the rims), because the GeV electrons emit¬ 
ting the radio are not affected b y losses. This is indeed what is observed in 
Cas A ( Vink and Laming . 20031) . In SN 1006 ( Long et all 120031 Fig. 4), it 
is true as well that the non-thermal X-ray emission (above 1.2 keV) is much 
sharper than the radio emission. The radio is never more peaked than the X- 
ray emission below 0.8 keV, which is mostly thermal (O K lines). In Tycho on 
the other hand, a radio pe ak just behind the blast wave (associated with the 
X-ray rim) seems to exist (Dickel et ah. 11991 ). This is difficult to understand 
in the synchrotron interpretation of the X-rays. 


Recently, Lvutikov and Pohli ( 20041) suggested that the narrow rims could 
mark the interface between the ejecta and the shocked ambient gas, rather 
than the blast wave. Indeed it is relatively natural to reach high magnetic 
fields there. However, I do not think this interpretation is tenable for three 
reasons: 

• Optical Ha observations show filaments at the same distance from the SNR 
center as the X-ray rims (although not coincident with them). The line 
shape indicates that a fast (originally ionized) and a slow (originally neutral) 
population exist there. It is hard to imagine what could cause this far from 
a shock. 

• No strong metal lines are observed immediately behind the rims, even in 
SN la (Tycho), where no hydrogen envelope should exist just behind the 
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interface. 

• One would have to find another explanation for the bulk of the radio emis¬ 
sion (far behind the X-ray rims), which is currently interpreted as originat¬ 
ing in the high magnetic fields at the interface. 


4 Physical conditions behind the blast wave 


Quantitatively, the interpretation of the X-ray rims as synchrotron goes as 
follows. The synchrotron cooling time depends on the magnetic field and the 
electron’s energy as t coo \ = 6.37 x 10 8 T>“q E ~ 3 s. Since the characteristic 
frequency at which an electron radiates is z/ sync = 1.82 x 10 15 B mG E^ Ig Hz, 
the cooling time may be expressed as a function of the frequency at which the 
rims are observed. 

f cool = 5.5 x 10 7 B~f s (1) 


Two effects (expected to be of the same order at the maximum energy reached 
by loss-limited electrons) combine to set the scale height of the emission be- 
hind the shock: advection jjBalletJ^ 200,'jj^ Vin k and Laming ._ 20031 Barnba et al. 
2003h and diffusion (|Berezhko et all 20031: Ya.may.aki et al 


Zadv = t coo iv sh /r = 1.8 x 10 3 bJ{ 2 i/ ke y /2 Uiooo/r pc 
kif = x/^dtcooi = 1-2 x 10~ 3 B~^ 2 pc 

_1 /o 

Ldv/^dif = 1-5 Z4eV Viooo/r 


( 2 ) 

( 3 ) 

( 4 ) 


where u s h is the shock speed, rqooo = u s h/(1000 km/s), r is the total compres¬ 
sion ratio and n d = c/ (3e) E/B is the downstream diffusion coefficient in the 
Bohm limit. 


The observed characteristics of the five young shell-type SNRs for which 
enough data exists are summarized in Table 1, together with the required 
magnetic field. The ratio / a dv/0if is always estimated assuming a compression 
ratio of 4. If the compression ratio is larger (as predicted by non-linear ac¬ 
celeration models), then Z ac j v will be smaller but the magnetic field estimate 
(from l di{ ) won’t change. 


Overall, this interpretation imposes B ~ 200 /iQ downstream in the three 
youngest SNRs. It is a very im portant result because it provides observational 
evidence for the idea (Bell and Lucek, 12001 ) that diffusively accelerated par¬ 
ticles streaming ahead of the shock are able to generate a turbulent magnetic 
field larger than the original ordered field (which cannot be larger than a 
few /j ,G in such surroundings). This is potentially the key for breaking the 
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Name 

Distance 

Shock speed 

Age 

Cut-off frequency 

Cas A 

3.4 kpc 1 

5200 krn/s 2 

320 yr ? 4 

1.2 keV 3 

Kepler 

4.8 kpc 5 

5400 krn/s 6 

400 yr 

0.9 keV 8 

Tycho 

2.3 kpc 9 

4600 km/s 10 

430 yr 

0.29 keV 11 

SN 1006 

2.2 kpc 12 

2900 km/s 13 

1000 yr 

3 keV 15 

G347.3—0.5 

1.3 kpc 16 

4000 km/s ? 17 

1620 yr ? 19 

2.6 keV 18 

Name 

Obs. freq. 

Zadv/Zdif 

Projected width 

Magnetic field 

Cas A 

5 keV 

0.88 

0.05 pc (3") 3 

230 n G 

Kepler 

5 keV 

0.92 

0.07 pc (3") 7 

180 ^G 

Tycho 

5 keV 

0.78 

0.05 pc (4") 11 

250 ^G 

SN 1006 

2 keV 

0.78 

0.2 pc (20") 14 

87 nG 

G347.3—0.5 

2 keV 

1.06 

0.25 pc (40") 18 

79 ^G 


Table 1 

Characteristics of the non-thermal emission behind the blast wave in young SNRs. 
The cut-off frequency is rather uncertain. The magnetic field is always estimated 
from equating Zaif (Eq. 3) and the projected width divided by 4.6 (App. A), even 
when Z a Hv/faif (from Eq. 4) is fo rm ally larger than 1. 


References: 1 iReed et all Jl99nh . 2 IVink et a ll d 19981). 3 Ivink an d Laming (l2003ll 


from a fit t o the spe c trum of the whole SNR, 4 lAshwort,hl ( 1 98(t l . 5 Revnoso and Gosd 
( 1999il. 6 _Hugh esl Jl999l)+d ista.nce. 7 Decourche lle (private c ommu nica tion), 8 


C^ggMmOignaLSi _al. _ l2004ali in the southeast, 9 ISmith et alJ dl9 9 ill. 10 

I Onnfjh 11 TTnrnn rr rvf oil JoflHO \ 12 ISTvnlrlov of ol ( 200yl ^ of- o 


Hwang et alJ (120021 

Bamba et ah 


Wi nkle r et al, 

( 200 : 


m 


Ghava mian et al. 

15 iRothenflug 


Hughei 

(^ooa) 


et al 


( 2f)nnh 

in the northwest, 

1 299411 in the northeast, 16 __ MM _ - ... , . 

ing the expansion nara,meter A = 2/3. 18 lBazenclic et ail 129931 in the northwest, 


in the nort heast. 

C assam- C heiuu_et_aL <200-1 1B. 17 from v = A r/t tak- 


19 


Wang et al, (1997) 


lLagage and Cesarskv (1983) limit and accelerating protons and heavier ions 
(not limited by radiative losses) np to the ’knee’ of the cosmic-ray distribution 
at 3 10 15 eV. An observational consequence is that it predicts weak inverse 
Compton emission in the TeV range, because the density of accelerated elec¬ 
trons (for a given synchrotron emission) is much lower than estimated from a 
’standard’ compressed held of 10 /iG or so. 


A somewhat puzzling aspect is the geometry of the acceleration with respect 
to the magnetic held. In G347.3—0.5 and Kepler, the very large asymmetry 
in the X-ray image is probably due to variations of the exterior conditions 
(ambient density), so it is difficult to infer anything on the geometry. In SN 
1006, the non-thermal emission is very clearly bipolar. The very faint X-ray 


emission from the center of the remnant above 2 keV led Rothenhug et al 


7 














































































(20Q4) to conclude that the non-thermal X- rays are e mi tted in two polar caps. 
This is coherent with theoretical ideas ( Volk et al. . 2003) that acceleration 
works better where the magnetic held was originally parallel (to the shock 
velocity). Unfortuna tely, that nic e picture is not at a ll what we see in Cas A 


( Hwang et al. . 2004 ) and Tycho ( Hwang et al. . 2002f) where the non-thermal 


X-ray rims seem to run all around the remnants. The main differences between 
those two and SN 1006 are the shock velocity and the ambient density (both 
lower in SN 1006). It may be that somehow the faster acceleration at larger 
shock velocities can work whatever the magnetic orientation, but a quantita¬ 
tive explanation does not exist yet. 


5 The shock precursor in the ambient gas 


The principle of diffusive acceleration predicts that the accelerated particles 
and the magnetic held should also be present some distance ahead of the 
shock. This means that synchrotron emission should be observed at some 
level upstream of the blast wave. 


The absence of such a precursor in the radio range (an upper limit on its 
width) implies an upper limit on the diffusion coefficient at the energy of 
electrons which radiate in the radio. This is equivalent to a lower limit on the 
turbulent held which is l arger t han the ordinary interstellar turbulence. This 
argument led Achterberg et al. (11994 ) to conclude that accelerated particles 
could indeed generate the turbulent held which is required for acceleration to 
be a fast process. 


Because the diffusion coefficient (for a given level of turbulence) increases with 
energy, X-ray observations (corresponding to electron energies 10 4 times larger 
than radio observations) provide a much more stringent constraint. In Cas A 
and Tycho, the radial prohle of the X-ray emission appears symmetric, so one 
might think that the width used above may be used as an upper limit to the 
size of the precursor. 


It is not so, though, because the shock compression necessarily results in 
compression of the magnetic held. Assuming isotropic magnetic turbulence 
upstream and shock compression by a factor r sub , the magnetic held down¬ 
stream may be larger than upstream by a factor re = ^/(l + 2rf ub )/3, and the 
synchrotron emission by approximately (because increasing B moves the 
spectrum both up and to larger frequencies), where T is the X-ray photon in¬ 
dex (typically 2.5 to 3). r sub here is the compression ratio at the gas subshock. 
It is always lower than 4, and decreases when the total compression ratio r 


increases. 
























In SN 1006, there is no evidence of a precursor either. Long et ah (2003) 
give a very stringent upper limit of 1.5% of the postshock level to the X- 
ray brightness about 20"(or 0.2 pc) upstream. They argue that this could 
be explained with a large magnetic field jump (l arger than 4j, p ossibly due 
to magnetic field amplification at the shock. Bamba et al.1 (2003) argue that 


it can be d u e to a perp endicular (to the shock normal) magnetic field. But 
Rothenflug et al. ( 20041 . on the basis of the overall geometry of the X-ray 


emission, show that this cannot be true (the bright limbs must be polar caps 
rather than an equatorial belt). iBerezhko et al. (2003! show that the observed 
profiles are compatible with the expected emissivity jump for a turbulent field 
(r B = 12.5 with r su b = 3.6 and T = 2.3), plus the intensity decrease outwards 
on the diffusive scale height. In other words, this constrains the upstream 
diffusive scale height for electrons emitting 1.5 keV X-rays to be lower than 
0.12 pc or so. Assuming Bohm diffusion, the upstream scale height ft u /v s h 
deduced from Table 1, dividing B by r B and deducing E corresponding to 
emission at 2 keV, is 0.08 pc. This is close enough that the precursor should 
be detectable in the near future. 
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A Projection of a thin sheet along the line of sight 


Let r be the distance to the center of a sphere in units of the sphere’s radius 
R s , and p the distance to the center of the disk it projects onto a plane (in the 
same units). The projection of a spherically symmetric emissivity S(r ) onto 
that plane results in an axisymmetric brightness profile 

\/i-p 2 

B(p)=2R s J S(r)dz 
o 

2 2,2 
r —p+z 

where z is the coordinate perpendicular to the plane. 

If the emissivity is concentrated near r = 1, it may be written £(r) = f(x) 
where x — (1 — r)/a, f is the functional form of the decrease and a is the 


(A.l) 

(A.2) 
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Fig. A.l. Emissivity profile exponentially decreasing inward (Eq. A.8 with a = 0.01, 
dashed line) compared with brightness profile projected onto the sky (Eq. A.l, solid 
line). The brightness units are R s x the emissivity at the shock, where R s is the 
shock radius. On that scale, the brightness at the center of the disk is 0.02. The 
emissivity profile was scaled down (for illustration purposes) to the same maximum. 
The dotted line (barely distinguishable from the solid one) is the res ult of the thin 
sheet approximation (Eq. A.5). The dash-dotted line is Eq. 6 of Berezhko and Volf] 

J23). 


scale height (in units of R s ). Defining y = (1 — p)/a and u = z/\/2a, r may 
be developed for p close to 1 at order 2 in powers of z to give 


r ^ p { l + jp) (A - 3) 

- ~y — u (A.4) 

a 

Vv 

B(p) ~ 2 R s V2a j f(y — u 2 )du (A.5) 

o 

Eq. A.5 is of the form B(p ) = 2R s y/2a g(y). It always has the same functional 
form in terms of y near the edge of the disk (whatever the width a). In the 
same limit of small a, the brightness toward the center of the sphere is 
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(A.6) 


1 

B(0) = 2R S J S(r)dr 
0 

oo 

£>(0)~2A s a J f(x)dx (A.7) 

o 


The formulae above may be quantified in the special case of exponential de¬ 
crease, shown on Fig. A.l. 


£{r) =exp 



1 


y/y 

B(p) ~ 2 R s V2a e~ y j exp (u 2 ^j d u 
o 


(A.8) 

(A.9) 


Then the maximum of g(y) occurs at y o = 0.854 where it reaches go = 0.541. 
g(y) decreases inwards to half that value at y\ = 4.685. g(y) decreases outwards 
to half that value at y 2 = 0.082. The full width at half maximum is thus 
(y i — 2 / 2)0 = 4.603 a. For comparison, the FWHM of the exponential itself 
is log 2 a = 0.693 a. The ratio between the brightness at the center of the 
sphere and at maximum is £>(0)/£> max = Ja/2/go = L307 y/a. The explicit 
(non-integral) approximation proposed by Berezhko and Volk (2004) is less 
accurate in representing the peak, but correctly represents the profile toward 
the center of the sphere (Eq. A.9 doesn’t) . Note that the difference with the 
value 7 a given for the width by Berezhko and Vo lk! ( 2004 ) is mostly due to a 
different definition (they consider the width at 1/e of the maximum). 
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